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Abstract: The dimensional synthesis is done by using Genetic 
Algorithm to achieve a desired trajectory. Three problems are 
analyzed having different curvature. The program is authored 
in MATLAB® 2010a. The error is seen to be in the permissible 
prescribed limit. The prototyping of straight line trajectory 
analysis is also done in ADAMS®. 
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1. Introduction 

The synthesis of mechanisms is usually employed to project or 
predict the functional dependence between geometry and the 
motions of various parts or links in mechanism under 
consideration. Conventionally the synthesis is classified into type 
synthesis, number synthesis and dimensional synthesis. The 
dimensional synthesis means determination of the dimensions of 
the individual links. Further the path generation is acquired by 
either with prescribed timing while the angle of each prescribed 
point of coupler curve are required to be worked out or without 
prescribed timing while the angle of each prescribed point of 
coupler curve is given. 

The classical graphical and analytical methods for the 
dimensional synthesis are discussed in Theory of Machines and 
Mechanisms by Joseph Edward Shigley(1988)which are 
restricted up to five precision points wherein the coupler curve is 
defined.In the case of more than five points the Genetic 
Algorithm (GA) is implemented for dimensional synthesis. The 
GA is an optimization technique which is based on natural 
evolution theory and follows the Darwin's principle of survival 
of fittest. Professor John Holland investigated the concept of GA 
in mid 60's for explaining the adaptation process of natural 
system for creating artificial system that works similarly. The GA 
converts natural evolution to world of computers, to solve 
optimization problem. Individuals those are weak will die off, 
and fit individuals will be selected for reproduction. 
Cabrera et al. (2002) formulated the GA and applied to optimize 
the Euclidean distance error function.This error function is the 
sum of squares of the difference between the desired path p~oint 
and generated path point. The intension is to reduce the error 
between desired point and generated point after every iteration 
so that at the end of last iteration the error should be within 
prescribed limits. Laribiet al. (2004) formulated fuzzy logic 
based GA and optimize the error function called as the 
orientation structural error of fixed link. Acharya et al. (2009) 
also implemented different evolutionary algorithmsvizGA, 
differential evolution and particle swarm optimizationfor the 
synthesis of four bar mechanism. The intension is to compare the 
efficiency of different algorithms. Cabrera et al. (2010) again 
implemented an algorithm which is the modification of the 



previous algorithm and the name of this algorithm is Malaga 
University Mechanism Synthesis Algorithm (MUMS A). 
Matekar et al. (2012) formulated the modified distance error 
function based on longitudinal and transverse errors between 
prescribed path points and obtained path points. 
In this paper the GA is formulated consisting of three operators 
namely reproduction, crossover and mutation and Euclidean 
distance error function is used for optimization. The algorithm is 
implemented on straight line generation problem passing 
through six desired points. The algorithm is implemented in such 
a way that, after application of every operator the population size 
is updated; so population size is updated three times in one 
iteration. 

2. Formulation of Objective Function 

The Fig.l shows four bar path generating mechanism O^ABO^. 
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The position of point P (Px., Py^) is calculated at each desired 
point i from above equation and then given as an input to the 
equation of error function (1) to evaluate the error. 
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Fig.lConstruction Details of Four bar Linkage 

2.1 Objective Function 

The objective of path synthesis problem is to have a minimum 
position errorand to achieve the objective, dimensions of the four 
bar linkage i.e. the design variables namely a,b,c,d,e,f links, 
input angle O2 and 6^ angle measured with the fixed link d are 
determined in such a way that the error between specified point 
{Pj^^. , Py^. ) and generated point (P^j , Pyi ) should be minimized, 

where i = 1,2,3, ,n and n is number of prescribed points. 

The optimization function is Euclidean distance error function is 
given below, 

fix) = Sill [ (P.,, - P^iY + {Py,, - Pyif] (1) 

where, X is the set of design variables to be obtained by 
optimizing error function. 
X=[a, b, c, d, e, f, ^1,^2]. 

2.2 Constraints 
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The significant equations in the synthesis of mechanism, (1) the 
objective function and (2) the equality or the inequality 
constraints which are applied to get optimum solution. The 
constraints include: Grashof s criterion, input link order, the 
transmission angle constraint, while the input link should be the 
smallest and the boundary value of variables. 
Grashof s criterion: 

This criterion is the prevailing condition for the mechanism to 
make the complete revolution of crank or the input link. 
s + I <p-\-q, where, 5 and / are shortest and longest links 
respectively and p and q are the other two links of the 
mechanism. 

Sequence of input angle constraint: 

A6/2^ = (62' - 6/2'"^), where i is the position of prescribed 
point. 

Transmission angle constraint: 

In order to get smooth link movements with respect to each 
other, the force and motion is required to be effectively 
transferred between them, which prevails that the transmission 
angle is close to 90°. This is a constraint for the designer, to 
maintain it in between allowable minimum and maximum 
values. 

Mmax = COS [ — J (2) 

Mmin = COS [ — J (3) 

The actual value of transmission angle at any angle 62^ is given 

by, 



/i = COS 



-1 



b^- - + +2adcos 



2bc 



(4) 



The criterion to be satisfied is /i^^j^ < /i < 
The length of Input link must be small: 

The length of input link must be smaller than the lengths of the 
remaining three links. This is strictly observed in the initial stage 
so that the constraint is not violated and the penalty is 
automatically written off. 
Variable bounds: 

Minimum and maximum values of all variables in the design 
vector are considered in the initial set of random variables and in 
order to get the value of variable between its minimum and 
maximum value, the simple generation rule is,X = X^i^ + 
rand^(X^ax -^mm)(5) 

where, rand is a random number generator between 0 and 1 . 
Overall Optimization Problem: 

The objective function is the sum of the error function and the 
penalties added on violation of the constraints as follows; 
E(X) = f{X) + PI * C1{X) + P2 * C2{X) + P3 * C3{X) 
(6) 

where, f{X) is the Euclidean distance error function. C1(X), 
C2{X), C3{X) are multiplying factors for Grashof s criterion, 
sequence of input angle and transmission angle criterion 
respectively and if the criterion is satisfied then corresponding 
multiplying factor is taken as 0 otherwise 1. PI, P 2, P 3 are the 
penalties which are to be added in the function value if these 
constraints are not satisfied and these penalties are taken as 1000 
as surveyed in various literatures. 

3. Optimization Algorithm 

The flow chart of applied algorithm is given below 
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step 1: Initialize Parameters 



Generations = No. of generations 
Popsize = Population size 
Cp = Crossover probability 

Mp = Mutation Probability 

L = Real value to control the 
crossover disturbance 

rangemu = Range for mutation 

X = Set of design variables, 



Step 2: Initialize X 



for i = 1 to Popsize do 



Ftandomlv generate the 
design variables 



Step 3: Application of constraints 



step 4: Reproduction 





fori = 1 to Popsize 




- Calculate /(JTJ 


i 



fori = ltoPopslze 



Evaluate the fitness value, F^^ 



Average fitness, Xr^ = suin(i^} /Popsize 



Expected Count 




for i = 1 to Popsize 




1 i 






Expected count, E;ri = f^ 







Probability of Selection 



for i = 1 to Popsize 



Update the matrix X as per the probability 



T 



. = Best individual vector 
^ri' ^rz = Randomly selected two individuals 




I Update Xj I 



Step 7: Termination Criteria 



f or i = 1 to Generations "| »|^^^p3j i|jtep4j *|jtep 5 



Fig.2Flow Chart of GA 
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4. Results and Discussion 

The key parameters of an algorithm are efficiency and accuracy. 
The problem statements given in this paper is to generate a 
straight line passing through 6 desired points without prescribed 
timing 

4.1 Straight Line Generation Problem 

Design variables are, 

X = [a, b, c, d, e, f,X^, Y^, 6^, ^2'. ^2' ^2'] 

where, a, b,c,d = Links of the mechanism 

e, f= Construction lines for obtaining the coordinates of 

coupler point P 
61 = Angle made by fixed link d with X-axis 

62 ^ = Angle made by input link with fixed 

link at i* position 
Target points chosen: 

C/ = [(20,20), (20,25), (20,30),(20,35), (20,40), (20,45)] 

Limits of the variables: 

a,b,c, de [0,60]; cJ^X^JqE [-60,60]; 

^1, 02', ^2'. ^2'. ^2^ ^2'. ^2'^ [0,27r]. 

Parameters considered are. Generations^ 1000, = 100, Cp = 
0.9, Mp = 0.1, L = 0.6, range = 0.1. 



Table. 1 Synthesized Results for Straight Line Generation 
Problem 



Design 
Variable 


Value 


Design 
Variable 


Value 


a 


10.4657 


62' 


5.33582 


b 


26.4562 




5.52386 


c 


43.3573 


62' 


5.69795 


d 


48.3378 


62' 


5.86718 


e 


-9.9783 


02 = 


6.03998 


f 


-47.2007 


02' 


6.22731 


^0 


57.6963 




0.09769 


Yo 


46.1450 







The synthesized values of all 15 design variablesinvolved are 
shown in Table. 1. The corresponding values of desired points 
(Px^, Py^) and generated points (P^, Py) and also the error 
between them is shown in Table. 2. 



Table. 2 Actual Points Traced by Coupler Point for Straight 
Line Generation Problem 



Sr. 
N 
0. 






-Pxf 


Pya 


Py 


^Pya 
-Pyf 


P.Y + 
(Pya 














-Pyf 


1 


20. 


19.96 


0.00138 


20. 


19.99 


1.2602 


1.40E- 


0 


281 


3 


0 


645 


E-05 


03 


2 


20. 


20.00 


7.34E- 


25. 


25.00 


4.761E- 


7.34E- 




0 


857 


05 


0 


007 


09 


05 




20. 


20.02 


0.00072 


30. 


29.99 


8.41E- 


7.22E- 


3 


0 


686 


2 


0 


997 


10 


04 




20. 


20.03 


0.00150 


35. 


34.99 


3.1541 


1.51E- 


4 


0 


88 


5 


0 


822 


E-06 


03 




20. 


20.02 


0.00089 


40. 


39.99 


6.561E- 


8.99E- 


5 


0 


997 


8 


0 


919 


07 


04 




20. 


19.93 


0.00424 


45. 


45.00 


1.4859 


4.25E- 


6 


0 


483 


8 


0 


122 


E-06 


03 
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The path traced by coupler point is shown below in Fig. 3. 
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Fig. 3 Path Traced by Coupler Point for Straight Line 
Generation Problem 
The objective is to generate a straight line trajectory from a 
coupler point of the four bar mechanism. The coupler point is 
the point on a triangular coupler link. The four bar mechanism 
contains four links namely input link, triangular coupler link, 
output link and fixed link as shown in Fig. 3. The coordinates 
(20,20) and (20,45) are the starting and end points of the desired 
straight line. There are 6 desired points through which the 
coupler point should pass. The generated coupler curve is shown 
in Fig. 3 and it is in elliptical shape but the right portion of 
generated elliptical curve is straight line and this straight line is 
passing through all desired points. The total error between 
desired and generated points is 0.0088. 



Generation Plot for Problem 1 
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Fig. 4 Plot of Error Vs Generations Count for Straight Line 
Generation Problem 



The GAimplemented in MATLAB®2010a gives off the graph of 
error verses generations count is plotted as shown in Fig. 4. 
Initially the structural error is more than 6000 and it is reduced 
up to 1 after 100* generation. 1000 generations carried out and 
at the end of 1000* generation the structural error reduces to 
0.0088. The error in X and Y direction at each position is 
calculated and plotted as shown in Fig. 5. 
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Fig. 5 Error in X and Y Direction for Straight Line Generation 

Problem 
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Simulation of mechanism: The prototyping of mechanism is 
done in ADAMS® to co-relate the simulation results with 
generated curve. The dimensions obtained fromthe synthesis of 
straight line analysis problem are used to form a mechanism in 
ADAMS®. The mechanism consists of input link, coupler link, 
and output link. The circular motion of input link is converted 
into elliptical motion of coupler point of triangular coupler link. 
If the input link is rotated in anticlockwise direction then in the 
return motion, coupler point gives straight line and it is 
analogous with the generated coupler curve in MATLAB®2010a. 
The coupler curve traced in ADAMS is shown in Fig. 6. 
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Fig. 6Coupler Curve traced in ADAMS 
5. Conclusion and Future Scope 



In straight line generation problem, the error is reduced to 0.008 
which is because population size is updated three times in each 
iteration. The obtained coupler curve is analogues with the 
desired curve. The OA allows more number of design variables 
as compared with analytical and graphical methods. 
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